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SUMMARY 

A finite-difference method 'previously presented for comput- 
ing elastic stresses in rotating disks is extended to include the 
computation of the disk stresses when plastic flow and creep are 
considered. A finite-difference method is employed to eliminate 
numerical integration and to permit nontechnical personnel to 
make the calculations with a minimum of engineering super- 
vision. Illustrative examples are included to facilitate explana- 
tion of the procedure by carrying out the computations on a 
typical gas-turbine disk through a complete running cycle. 

The results of the numerical examples presented indicate that 
plastic flow markedly alters the elastic-stress distribution. 

INTRODUCTION 

With the advent of jet propulsion as a motive force for air- 
craft, the gas turbine has become an important source of 
power. In most machinery, design stresses are limited by 
the yield strength or the creep strength of the material em- 
ployed, together with a certain factor of safety, and little or 
no analytical consideration is given to the occurrence of 
plastic flow under operating conditions. Gas-turbine disks, 
however, are required to operate under thermal gradients 
and centrifugal forces producing stresses that, in materials 
currently available, frequently exceed the yield strength, 
resulting in plastic flow. The interaction of plastic flow and 
creep, together with the variation of thermal gradients 
through a series of cycles consisting in starting, running, and 
stopping, can produce stress distributions and even failures 
that might not be anticipated on a basis of elastic-stress 
analysis. 

A rapid routing method of elastic-stress analysis of rotating 
disks is presented in reference 1 , which gives accurate values 
of the true stresses in disks provided that the yield strength 
of the material is not exceeded. The finite-difference method 
of reference 1 has been extended at the NACA Cleveland 
laboratory to include consideration of plastic flow and creep, 
which thus allows calculation of the true stresses in a gas- 
turbine disk and gives the variation of stress distribution 
with time. The handling of plastic flow is somewhat less 
routine than the calculation of the elastic stresses in that a 
repetitive trial procedure is required. With practice, the 
correct value can be obtained on the fourth or fifth trial. 
The computation of the effect of creep, although in procedure 
the same as the computation of plastic flow, is a direct calcu- 


lation requiring no trial-and-error procedures. Because the 
method eliminates numerical integration, nontechnical per- 
sonnel can make the calculations with a minimum of engineer- 
ing supervision. 

SYMBOLS 

The following symbols are used: 
c creep rate under stress a„ (in./(in.)(hr)) 

E elastic modulus of disk material, (lb/sq in.) 
h axial tliiekness of disk, (in.) 

R rati0 (fe) 

r radial distance, (in.) 

T temperature, (°F) 

u radial displacement, (in.) 

a coefficient of thermal expansion between actual tem- 
perature and temperature at zero thermal stress, 
(in./(in.)(°F)) 

T total creep under stress e>, (in./fin.) 

A plastic increment of strain, (in./in.) 

A r plastic increment of strain in radial direction, (in./in.) 
A, plastic increment of strain in tangential direction, 
(in./in.) 

A T temperature increment above temperature of zero 
thermal stress, (°F) 

S r creep increment in radial direction, (in./in.) 

S t creep increment in tangential direction, (in./in.) 

e strain, (in./In.) 

e, plastic strain corresponding to stress <r« in tensile speci- 
men, (in./in.) 

€ r radial strain, (in./in.) 

et tangential strain, (in./in.) 

fi Poisson’s ratio 

p mass density of disk material, ((lb) (sec 2 ) /in. 4 ) 
cr stress, (lb/sq in.) 

<r t equivalent tensile stress, (lb/sq in.) 

£r r radial stress, (lb/sq in.) 

o', tangential stress, (lb/sq in.) 

<r v proportional elastic limit, (lb/sq in.) 
t time during which creep occurs, (hr) 

a angular velocity, (radians/sec) 

The following supplementary subscripts are used for de- 
noting values of the preceding symbols in connection with 
the finite-difference solution: 
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n nth point station 

7i-l (n-l)st point station 

a station at smallest disit radius considered 

(For disk with a central hole, this station is taken 
at the radius of the central hole; for a solid disk, this 
station is taken at a radius approximately 5 percent 
of the rim radius.) 

b station at rim of disk or base of blades 


The following supplementary symbols denote combinations 
of the foregoing symbols: 

A r , a 

A l.n 


-Br.it 
Bt,nj 
C a =r„h n 


stress coefficients defined by equations 

}'0'r,n — A r<n <T + Pr,n 

& (,tt = A [ tn (T Pi.it 


m flit i (1 -|~ ^n) 0\< l) 

° 2 E n r n 

En (r n r 


D\ 


1 , (1 + Mn) (r n —r n _i) 

= F S + 2E7 n 


F n — T n -ihn-\ 


pi/ fln-l (1 T Mn-l) 0"n ? i»-l) 

’* Eh- i 2E n -ir a -i 

G n =2 0'n 1'n— l)b n —\ 

pr 1 (t ~H Vn-l) (?*» ?n— l) 

n E n ~i 2En-ir n -i 

H n — g W ! (r B r„_i) {pnhnTnA- Pn-^-n-^n- 1 ) 
H' n =a n AT„—aH-iAT n - 1 


„ F'JDh-FhD'h 
K i! ~C\D n -C n D\ 

K \ ?ifL n - c ' nF , n 

C' iiDn CJDh 

r G'hDh+GhD'h 
- L " ~ C'hDh-ChD'h 


T / C'hGh+ChG'h 
V'hDh-CJ)^ 



,, D'hHh+DAH'h-P'h-Q'h) 

M »- TFJJ^UJK 

w, C'JI n +C n (H' n -P\-Q ' n ) 

m " CnDn-CnD’H 

(A/„ and M‘ n are defined in reference 1 for the special 
case P' n =Q' n =0 ) 


O-’saO 

+ ,j! i h+ 


Sf.n-i 


0 


r»— ?V-A 

2r„_! / 


ANALYSIS OF PLASTIC FLOW AND CREEP 

Assumptions. — Four assumptions are made in tho subse- 
quent analysis: 

1. The disk material is linearly elastic up to a limiting 
stress value, called the proportional elastic limit, and above 
this limit plastic flow occurs. 

2. All variables of material properties and operating 
conditions are symmetrical about the axis of rotation. 

3. Axial stresses may be neglected and tho radial and 
tangential stresses are uniform across the thickness of the 
disk. 

4. Temperatures are uniform across the thickness of the 
disk. 

Outline of method. — In any thin rotating disk, the com- 
plete stress state is defined when the two principal stresses, 
radial o> and tangential <r ( , are known at every radius. Two 
equations relating these stresses to the radius arc required 
to specify the stress distribution. The first of these equations 
can be determined from the conditions of equilibrium of an 
element of the disk and involves no elastic properties of the 
material. The second is derived from the compatibility 
conditions, which state the interrelation of radial and 
tangential strains. The compatibility conditions arc depend- 
ent upon stress-strain phenomena and must thcreforo 
includo any departure from linear elasticity. When modifica- 
tion to allow for any possible departure from Hooke’s law 
is made, the compatibility conditions become true for any 
value of stress. The equation derived from the compatibility 
conditions thus modified, together with the equilibrium 
equation, is treated by the finite-difference method of refer- 
ence 1, and similar equations are obtained. These equations 
result in additional terms in the final equations, which are 
used to modify the result of the elastic calculation. 

Whenever stresses under discussion have been calculated 
by the method of reference 1 only, they will be referred to 
as “elastic stresses”; where plastic flow and creep have been 
taken into account, the stresses will be referred to as “plastic 
stresses.”. 

Derivation of method. — Tho equilibrium equation, which 
applies to both the elastic and plastic conditions, is 

^(rho- r )—h<r t -\-p<i> i r i h= 0 (1) 


The elastic compatibility relations given in terms of the 
radial displacement are 


and 


dU <J r — 1*17 1 I 


a A T 


( 2 ) 



<T t — MO> 

E ' 


{■aAT 


(3) 



DETERMINATION OF STRESSES IN GAS-TURBINE DISKS SUBJECTED TO PLASTIC FLOW AND CREEP 


Equations (2) and (3) must be modified to include consider- 
ation of plastic flow. When a material is stressed beyond the 
proportional elastic limit, the strain in the material is 
different from that indicated by Hooke’s law. The strain 
under such a load may be considered as being made up of 
two components, one elastic as predicted by the laws of elas- 
ticity and one an increment of strain due to the flow that 
occurs. Rewriting equations (2) and (3) on this basis gives 


du 

* r= lr~ 

-_^J^ +aA T +Ar 

(4) 

u 

CT, IJ-Or . . rp 1 . 

/E\ 

6l ~r~ 


\y) 


Similarly, any creep that occurs represents an additional 
departure from elastic behavior, which further modifies 
equations (2) and (3) to 

( 6 ) 

*<=—=■ ‘ jg r +aA!T+A t -|-Si (7) 

When the parameter u is eliminated as in reference 1, 

( 8 ) 

Applying the finite-difference method to equations (1) 
and (8) and using the notation introduced in the section 
Symbols gives 

CVr.B — &n<r t,a=F n (r (?»<r ZZ» (9) 


The solution of the equations is facilitated by the substitu- 
tion of the stress coefficients .A,*, B T , n , and B t , n into 
equations (9) and (10) . Proceeding as in reference 1 results in 
the equations 

C n A rra — D n A t . a — F % A r , a -i — (?».d.j,*_i=0 

C' n A r , a -D f a A lM -F' 1l A T . a . 1 +G' K A t , a - 1 =0 

C\B r , a -D' n B t , ll -F' a B r , a . l + 

G'nB lra -i—H'n+ P' K -\- Q / » = 0 - 

All but the last of equations (11) and equations (15) of 
reference 1 are identical. When equations (11) are solved 
for d rn , At,n f Br w n] and Bt, a > 


A Ti & — isT,n r .,_ 1 + 

Ai iTt = K n A Tt1t _i-{- B aA^n—l 

B r ,a=K n B rtn _ 1 -\-B % B ttn ^i-\-M n 

B Ut =K?JBr.,- l +L' m B u >- l +M'. 


The symbols K n , K' n , L„, and L' n have the same meaning- 
as in reference 1 . The 21/* and M' n terms are now defined 
as 

D'nH n +D a (H' a -P\-Q\) 

C\D a -Cj}\ (13 > 

_ C\H n +C n (H' a -P\-Q\) ~ 

M n ~ C\D a -C tt D\ (13a> 

The elastic case of reference 1 thus becomes a special case of 

the more general problem in which P'» and Q' n are both zero. 

Evaluation of plastic terms. —In order to apply the finite- 
difference method to problems involving plastic flow, 
a relation between stresses and strains in the plastic region 
must be established. In reference 2, a numerical-integration 
method for computing disk stresses is presented in which 
elongation is assumed to proceed at constant stress when the 
proportional limit is reached. References 3 and 4 present 
equations for the plastic relation of stress to strain based on 
the maximum distortion theory. Relations can be derived 
from the equations given in reference 4, which form a con- 
venient means of finding the plastic increments corresponding 
to the stresses present in the disk. Rewriting these equations 
in the notation of this report and letting the subscripts- 
1, 2, and 3 denote the three principal directions in the most 
general case give 


A i — 3- [(0-1 — o-j) + (ot— tr 3 )] 

d 2 = ^ [(o-j — 0 -,) + (oj 0 - 3 ) ] 

J? 

Aj=g- [(ctj — °"i) d - (o"a — 0a)] 


v e =-L V(or— trj) J + (cr a — <r,) a -f (a*— cr x ' 


where the ratio i? is defined in terms of the corresponding 
uniaxial stress a e and plastic strain e„ in a tensile specimen 
by the relation 

(16) 

When equations (14), (15), and (16) are reduced to the biaxial 
condition, which is assumed to prevail in the disk (that is, 
cr 3 =0), and the finite-difference notation is introduced 


A/.* — g- («T r .» — O ’ (,*”{- 0V.it) 

Aj,*=g- (ca* — 

o’ «.»= Vv r.a 1 — O’ r .»ff t,»+ <T 


73 
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Figure 1. — Unlaxfal stress-strain carve showing relation between stress tr,,« and plasflo 

strain 


Substituting equation (19) into equations (17) gives 


\,n=K~ (2a> >n — <r (,„) 

A,„=^ (2ff t ,n — (Tr.re) J 


( 20 ) 


A typical uniaxial stress-strain curve illustrating the rela- 
tion between effective stress <r tM and effective plastic strain 
e„, R on such a curve is shown in figure 1. Investigations at 
the National Physical Laboratory of Great Britain on 
turbine-disk alloys and experiments by Taylor and Quinney 
(reference 5) were found to correlate well with equations 
(20) for the range of strains over which the volume of the 
material is approximately constant. 

Equations (20) give the relations between plastic strains 
and true stresses that will be used as the basis for numerical 
calculations in the present report. The method of stress 
analysis to be presented does not depend, however, on the 
validity of these equations. As more accurate relations are 
determined between stresses and strains, these relations may 
readily be used in place of equations (20). 

Calculation of plastic flow when no previous plastic flow 
has occurred, — The determination of the plastic stresses in. 
the disk resolves itself into the problem of finding corre- 


sponding stresses and strains that satisfy equilibrium and 
compatibility equations (9) and (10), and biaxial stress-strain 
equations {20). The problem is approached by first com- 
puting the elastic stresses, and the equivalent uniaxial tensile 
stress at each station is determined from equation (18). 
If at any station this stress exceeds the proportional clastic 
limit of the material at the temperature at this station, then 
plastic floiv takes place, and it becomes necessary to resort 
to a trial-and-error procedure to adjust the stresses to allow 
for this flow. 

Assume, for example, the equivalent uniaxial stress at a 
given station lies at point A on the extension of the modulus 
line in figure 2. Because the point A lies above the propor- 
tional elastic limit (point B), plastic flow must occur. The 
stress and the strain must he adjusted to fall on the curved 
stress-strain curve that is characteristic of the material. 
As a starting pointythe total strain in the true stress-strain 
condition is assumed equal to the strain at A. The stress- 
strain condition at the given station then lies on the constant- 
strain line through B, or at~C. The plastic strain is 



Fiouke 2. — Uniaxial stress-strain curve Illustrating procedure used to And correct value of 

plastic strain. 
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given by CD. Values of A r> „, A,,,, and P' n may be obtained by 
using this value of together with the values of ov,*, 
and <t, iK from the elastic calculations. 

Once P' n has been calculated, new values of cr rf „ <r t .n, 
and <r» tS can be computed. The new value of is greater 
than that at point D, such as that at point E. Although the 
stresses corresponding to tr e a at point E together with the 
strain CD meet the conditions of equations (9) and (10), 
they locate the stress-strain point F, which is not on the 
stress-strain curve, so that the physical conditions imposed 
by the material are as yet unsatisfied. Inasmuch as any 
value of t Pf n less than CD would give a value of a e .n greater 
than that at E, CD is a lower limit of Similarly, 

because the value of a, tK calculated by using an of CD 
is too great, the increment of strain EG corresponding to 
this <r e , a is an upper limit of Inasmuch as the true 

value of lies between CD and EG, their numerical 
average, shown as HK, is assumed to be a good approxima- 
tion. New values of P'n, <r r n , <r, n , and <r <n can be com- 
puted by using HK for €p n , the stress at E for a- en , and 
c T X and ar t ,n- Assume that this new value of tr 4 „ lies at 
the point M. Because the stress at M is higher than the 
stress at H in value, the increment H K is too small a value 
of e Pr „ and is therefore established as a new lower limit of 
6 p x . Further, because M is less than E, the corresponding 
increment MN is a new upper limit for e px and the process 
could be repeated again with the numerical average of M N 
and HK. Similarly, if the calculation using an t p „ of HK 
had resulted in a a en at P, HK would constitute a new 
upper limit and PQ a new lower limit. Had the resulting 
c e ,m been at R, HK would have still become the new upper 
limit of e Pi „, but CD would have remained as the lower 
limit. The process is repeated until the value of used 
in the computation and the e Pi „ corresponding to the 
resulting <r >iB are equal. 

Calculation of plastic flow when previous plastic flow has 
occurred. — The equations for strain that would apply to a 
disk that had already undergone the plastic strain are 

e,= ^y- £ +«AT+[A r ] + A r (21) 

s < = T^ r +“A? 7 +[A,3+A £ (21a) 

Here the terms [A r ] and [A,] represent strains already 
existent in the material before the application of stresses o- t 
and o> and are constant for the calculation, whereas A, 
and A t represent the components of plastic strain resulting 
from the application of c> and <r £ . In the solution of the 
equations by the finite-difference method, a term [P'„] 
appears together with term P'„. 'When previous plastic 
flow has occurred only once, [P'„] is identical with P' n from 
the previous calculation; where plastic flow has previously 
occurred more than once, [P'J is the algebraic sum of all 
earlier P' K terms. Thus, the previous plastic flow given by 
[P' n ] may be grouped with the temperature-effect term H' n 
by replacing H'„ with H' n —[P' J. 

This procedure amounts to an assumption that, as the 

003SSo— 50 10 



Figubk 3.— Uniaxial stress-strain curves showing components of strain when plastic flow 
occurs a second time. 

load and the temperature change, the stress position on the 
new stress-strain curve would be the same as if a test speci- 
men were loaded above the yield point, the load removed, 
the temperature changed, and a new load applied. This 
assumption is illustrated by figure 3, in which point A 
represents a loading at the first temperature condition; the 
dotted line AB represents the load-removal path; the curve 
BCD, the stress-strain curve at the new temperature; and 
point C, the new stress position. The total strain at point C 
is given by the sum of three strains. The residual strain 
caused by the first loading is ej is the elastic part of the 
strain caused by the second loading; and q, the plastic strain 
caused by the second loading. 

When the foregoing procedure is applied, the curve BCD_ 
must, of course, represent the true stress-strain curve at thfiL 
new temperature of a material that has already been sub- 
jected to the plastic cycle OAB. In general, the new stress- 
strain curve is different from the stress-strain curve at the 
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Figure 4.— Typical deformation-time curve from a constant-temperature, constant-load 

creep test. 

given temperature of a material that has not been subjected 
to plastic flow; however, unless data are available it may be 
necessary to assume that the curve BCD is the stress-strain 
curve at the given temperature of a specimen of virgin 
material. 

Calculation of effect of creep. — Creep is usually defined as 
the continuous deformation of material under a continuously 
applied load. Experimental data on creep of various ma- 
terials are usually obtained from tests run under constant 
load and temperature, although in many engineering appli- 
cations of materials the more general problem of changing 
load and temperature must be considered. The deformation 
curve obtained in a typical test is shown in figure 4. From 
this figure it can be seen that the deformation may be con- 
sidered as having occurred in three stages. During the 
primary stage, the deformation proceeds at a decreasing rate; 
during the secondary stage, at a constant rate; and during 
the tertiary stage, at an increasing rate, which proceeds until 
failure occurs. 

Because of the lack of data on creep except for uniaxial 
tensile stress, a relation between creep deformation and stress 
must be assumed. The following equations have been used 
for calculations in this report but, as better date, become 


Point stat ions 

a345 6 7 8 9 IQ U !2 13 14 15 16 !8 b 



.Fiauitz {.—Cross section of disk used for numerical examples showing location of point 

stations. 


available, more accurate relations can bo used. By the use 
of reasoning similar to that employed in determining the 
biaxial components of plastic-strain formulas for the creep 
increments, and 5, t » may be written 


««.»= 97 E -( 2 v t> »-v r . a ) (22a) 

*ve,n 

In equations (22) and (22a), r„ represents the total creep 
that would occur in time t under the uniaxial stress er«, R . It 
is here assumed that for sufficiently small values of r tho 
creep may be considered as occurring instantaneously at the 
end of the time period. 

During the secondary stage of creep, a characteristic creep 
rate c a exists,* corresponding to tho stress cr«,« at temperature 
T, and P s is given directly by 

r„=c»r (23) 

This rate is the value usually published in papers on creep 
and is the rate used for the numerical calculations of tliis 
report. "During primary and tertiary creep stages, the creep 
rate is also a function of time, but does not otherwise compli- 
cate the computation. 

Once values of 8 r jt and S t B have been found, the values of 
the Q' n terms may be determined and new values of cr r> « and 
may be computed. If the computed values of o>.» and 
<r,,» differ by more than a small amount, perhaps 2 percent, 
from the values of these stresses before creep occurred, 
a shorter_time interval should be selected and additional 
computations made for each such time interval required to 
equal the total time during which creep occurs. The effect 
of creep that occurred at previous time intervals is considered 
in a manner similar to that employed in considering previous 
plastic flow. The successive values of Q' n are summed to 
form a term [Q' H ], which gives the total effect of all previous 
creep deformation so that the term 

is replaced by 

H\-[P'n\-[Q'n] 

In any calculation of stress distribution subsequent to the 
occurrence of creep, the creep term [Q'„] is combined with 
the term [P'*] as the cumulative effect of all previous 
plastic deformation. 

Examples showing in detail how successive stages of 
plastic flow and creep are computed, each stage considering 
all previous plastic deformation, are given in the following 
section. 

NUMERICAL EXAMPLES 

The numerical examples presented hero represent a set of 
computations during one complete start-run-stop cycle for a 
typical turbine disk with a continuous rim and welded blades. 
The assumed profile of _the disk is shown in figure 5, together 
with the locations of the point stations used in the computa- 
tions. The assumed temperature distributions and corre- 
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sponding turbine rotative speeds are shown in figure 6. 
Curve IV and the corresponding speed of 11,500 rpm 
represent the steady-state running condition. Curves I to 
III and the corresponding speeds represent running condi- 
tions through which the turbine disk passes in reaching 
steady-state operation. Curves V to VII together with the 
respective speeds represent running conditions through which 
the turbine disk passes when being stopped. Creep is 
assumed to occur only during the steady-state running period. 

The physical properties of the disk material, including 
specific gravity, modulus of elasticity, stress-strain character- 
istics, and thermal coefficients of expansion, were based on 
the data appearing in reference 6, together with unpublished 
data obtained from the author of this reference. The stress- 
strain curves were constructed on the basis of these data and 
those for example I appear in figure 7. Inasmuch as no data 
were available on the effect of previous plastic flow on the 
shape of the stress-strain curves, it was necessary to ignore 
such effects and to use curves obtained directly from simple 
tensile-test data. Creep properties corresponding to a mate- 
rial having good creep resistance were assumed. The effect 
of primary creep was omitted because of lack of data. 

Because the disk used for these calculations is solid at the 
center, a supplementary numerical example showing the 
computation of the plastic-flow effect on stress distribution 
in a disk containing a central hole is given in the appendix. 

Example I. — Example I is the calculation of the stress dis- 
tribution in a disk operating under the conditions of curve I 
of figure 6 and having been subjected to no previous plastic 
deformation. These conditions are assumed to represent 
disk operation after the first short period of steady combus- 
tion when gas temperatures are high, thereby establishing a 
steep temperature gradient between the center and the rim 
of the disk. 

The preliminary elastic calculation is carried out in table 
I (a) by the method of reference 1. Two changes are made in 
the tabular setup. The first change is the insertion of 
columns 25a and 25b immediately following column 25. 
Column 25a lists the accumulated values of [P'*] and [Q'„], 
the total effect of previous plastic deformation. For the 
present example, this column is zero for all stations. Col- 
umn 25b is the value of the term H\— [P'J— which in 



Fig use 6.— Assumed temperature -distribution curves and corresponding turbine speeds. 
(Curves I, II, and in are consecutive starting conditions; carve IV represents steady-state 
operation; and curves V, VI, and VII are consecutive stopping conditions.) 


this example is the same as column 25. The second change 
is the computation of M K and (columns 31 and 32, 
respectively), which were computed in reference 1 by the use 
of column 25. In the present computations, column 25b 
is used. In addition, two more columns, 40 and 41, are 
added. Column 40 lists the values of the proportional 
elastic limit cr B . n of the material and column 41 lists the 
values of <r #>a as computed from equation (18). The entries 
in columns 40 and 41 of table I (a) show that the equivalent 
stress a-,'* is less than tr,,* for all point stations except 17 
to b. The effect of plastic flow must be considered at these 
stations and flow at these stations modifies the stresses at 
other locations in the disk. With the exceptions and the 
additions noted, the method of computation is the same as 
the method of reference 1 and wall not be discussed in further 
detail. 

The plastic-flow calculation has been divided into two 
parts because several quantities used in the computation 
depend only on the dimensions of the disk and can be used 
in all subsequent calculations involving plastic deformation. 
These quantities are computed for stations 17 to b, as shown 
by the four column headings of table I (b). 

The second part of the plastic-flow calculation is given in 
table I (c) . The first column in t his part of the table (column 
46) lists the values of e, a obtained from the corresponding 
stress-strain curve (fig. 7), as previously explained. Column 
46a lists the value of e p » used for the ensuing calculation, 
which for the first approximation is the same as column 46. 



(a) 4S0® F. (b) 623° F. 

fe) 809° F. CdJ 1050° F. 


Figure 7.— Stress-strain curves of disk material for various temperatures-. 
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Columns 47 and 48 list tlie values of A r ,„ and A/,*, respectively, 
computed by equations (20). Columns 49 to 52 are com- 
puted as shown by the column headings and from these 
columns the values of P r n are computed and listed in 
column 53 . . Column 54 gives the values of the term 
H'„~[P' „]-[(?' n]-P'n, which is then used to compute new 
values of M nt M\, B T Jl} B t n , a t ,a, ar,*, <r t ,n, and <r e , n as shown 
in columns 55 to 62, respectively. The new values of cor- 
responding to the new values of <r 4 ,» are read from figure 7 
and listed in column 46 of the second-approximation calcu- 
lation. The values in column 46 for the first and second 
approximations now constitute the lower and upper limits, 
respectively, of the possible strain increments. For the 
second approximation, column 4@a therefore lists as the 
values of e p ,„ to be used in this set of calculations the numeri- 
cal averages of the two sets of readings from the stress- 
strain curve. From this value, another new set of stress 
values is computed and a third set. of readings listed in 
column 46. 

At this point in the calculation, two alternate procedures 
are possible, as shown by consideration of station b. Inas- 
much as the average value of 4300X10 -8 inches per inch 
used in the second approximation gave a graph reading of 
1900 X10 -6 inches per inch, the averaging procedure would 

indicate that the next trial should.be X10~ 6 

or 4130X10 -6 inches per inch. This value could be used 
and the procedure continued until the. correct value is found. 
Considerable time may be saved, however, in making the 
calculation if a weighted approximation is used. Because 
the plastic-strain value of 3960 X10~ 9 inches per inch gave a 
resulting reading of 4650X10 -8 inches per inch whereas the 
value 4300X 10 -8 inches per inch gave the reading 1900X10 -6 
inches per inch, the strain 3960 X1Q7 9 inches per inch is 
apparently more nearly correct than 4300 X10 -8 inches per 
inch. In addition, the shape of the stress-strain curve in the 
region of 3960 X 10r 8 is such that smalLincreases in stress cor- 
respond to large changes in strain. If a trial calculation were 
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made using a value closer to 3960 than 4130 (for instance, 
4000), more information might be obtained than would be 
obtained by the averaging procedure. The right answer is 
thereby more quickly obtained. The same reasoning might 
be applied to the selection of values to bo used at the other 
stations for the third calculation. The second of these 
two procedures has been used in tablo I (c), as can be seen 
from the values of in column 46a used for the third 
approximation. 

Completion of the third approximation and comparison 
with new values of strain obtained from the stross-strain curve 
show the estimates of the third approximation to bo nearly 
correct, so that small adjustments made to compute the fourth 
and fifth approximations give the final answers. A calcula- 
tion equivalent to a sixth approximation is then made to 
column 53 to get the final correct values of the P' n terms. 
The stresses at the other stations a through 16 can now be 
computed by using the value of a t-a found in the sixth approxi- 
mation with the values ri r ,„, B r and B ti n found in 
table I (a). The values of plastic stress at all radii together 
with the elastic-stress distribution arc plotted in figure 8. 

Example n. — Example II considers tho disk studied in 
example I at tire time that the operating conditions have 
reached those indicated by curve II of figure 6, The elastic 
calculations are made by the method of reference 1 modified 
in accordance with the changes made in example I. The 
essential, parts of the computation are shown in tublc II, 
which is abridged from the complete calculation. Column 25a 
lists the values of [P' n ] that W'ere found as the final values 
of P\ in example I. Plastic flow occurs at stations 17, 18, 
19, and"6, and calculated true stresses when this plastic flow 
is considered are listed in table II. The stresses obtained 
as a result of this computation are plotted in figure 9. The 
elastic stresses obtained without considering the plastic flow' 
that occurred previously are also plotted for comparison in 
figure 9. 

Example HI. — Example III continues tho cyclo analyzed 
in examples I and II, at the conditions of curve III of figure 6. 


Figube 8.— Stresses in turbine disk under conditions of curve I of figure 0. 
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Figure 10. — Stresses In turbine disk under conditions of curve HI of figure 6. 


Table III gives the essential parts of the calculation for this 
example, which is similar to the procedure for example I in 
table I. The value of [P' a ], column 25a, however, is the 
total of the values of P'* obtained from examples I 
and II. In this example, plastic flow occurs only at stations 
17 and 18. The results of this computation, together with 
the elastic-stress curves found without consideration of 
previous plastic flow, are shown in figure 10. 

Example IV. — The steady-state operating conditions repre- 
sented by curve IV of figure 6 are treated in example IV. 
The essential calculations shown in table IV (a) were made 
similarly to those in table III except that [P'J in column 25a 
is now the sum of the values of P' a from examples I, II, and 
III. Because no values of <r,,» exceed those of a>,„, no plastic 
flow occurs and the stress values of table IV (a) are the true 
stresses at the beginning of steady-state operation. However, 
as parts of the disk are at elevated temperature, significant 
creep can occur at steady load at stations 16, 17, and 18 where 
the stresses are sufficiently high. Table TV (b) shows the 
calculations of creep. Column 63 lists the creep rate c n 
(in./ (in.) (hr)), and column 64 the creep increment r» for the 
5-hour running period. Columns 65 and 66 give the com- 
puted values of 5 r .« and respectively. The computation 




Figure 12.— Stresses in tnrbine disk under conditions of curve V of figure 5. 


then proceeds in a manner similar to the plastic-flow calcu- 
lations, as indicated in the column headings. The values 
for stress obtained indicate that, for small values, creep has 
only a slight effect on the stresses. Figure 11 shows the 
stress distributions at the beginning and the end of the 
steady-state running period, together with the elastic stresses 
obtained without considering either creep or previous plastic 
flow. 

Example V. — The conditions of example V represent one 
of the conditions through which the turbine disk is assumed 
to pass during the stopping period. The abridged elastic 
calculations are given in table V. Values listed now repre- 
sent the accumulated effect of plastic flow [P'„] plus the addi- 
tional effect of the creep represented by [<?'*] ; [$'*] is the same 
as the Q' n computed in example IV. All values of <r,,„ are 
less than the corresponding values of a-,,*; therefore no 
plastic flow occurs. The results of the calculation _ are 
plotted in figure 12, together with the elastic stresses com- 
puted without considering previous plastic flow or creep. 

Example VL — Example VI is the computation of the stress 
distribution at the temperature distribution assumed to be 
present shortly after the wheel has stopped turning. The 
essential parts of the calculation are shown in table VI. 
Because no flow was found in example V, the [P'J + IQ'J 
term will be the same in this example as in example V. 
Plastic flow occurs at station b. The resulting stresses are 
plotted in figure 13, together with the elastic stresses com- 



Tiotm U.— Stresses in turbine disk under conditions of curve IV of figure 6. (The stresses 
before and after creep occurs coincide within the accuracy of this plot.) 


Figure 13. — Stresses In tnrbine disk under conditions of curve VI of figure A 
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puted without consideration of previous plastic flow and 
creep. 

Example VII. — Example VII is the computation of the 
stress distribution in the disk after the temperature has 
bocorne uniform at the ambient temperature (assumed to be 
70° F) throughout, the disk. These stresses are therefore 
the residual stresses in the disk resulting from the flow 
occurring during the complete operating cycle. The abridged 
calculations given in table VII indicate that plastic flow 
occurs at station b and the residual stresses are plotted in 
figure 14. 

Discussion of numerical examples. — The foregoing cycle 
of stress calculations is indicative of the means of obtaining 
a complete analysis of the stress behavior of a turbine disk. 
Although the results plotted in the various figures and 
summarized in figure 15 do not represent the exact behavior 
of any particular turbine disk because of the lack of data on 
the material properties and temperature gradients, they do 
give a qualitative picture of. the behavior of a turbine disk 
with welded blades. The high residual tensile stress at the 
rim of the wheel provides a plausible explanation of the rim 
cracking that has occurred in such wheels. The compressive 
flow at the rim during starting and the tensile flow on stopping 
result in cyclic flow of the rim material with each start and 
stop and possibly induce cracks. When accurate data are 
available on creep, stress-strain relations, the effect of 
strain-hardening, and temperature distribution, quantitative 




Figche 15.— Plastic tangential stresses In turbine disk during one running cycle. 

analyses of disk behavior will be available as a guide in 
future turbine design. 

CONCLUSIONS 

A method for studying the operating stresses in gas-turbine 
disks has been presented that includes consideration of the 
effect of plastic flow and creep on the stress distribution. 
[Results of calculations indicate that rim cracking in turbine 
wheels with welded blade attachments may be caused by 
alternate compressive and tensile plastic flow as the wheel is 
alternately heated and cooled. From the results of tho nu- 
.merical examples presented, it may be concluded that plastic 
flow markedly alters the elastic-stress distribution. 

Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 

Cleveland, Ohio, March 5, 1943. 


Fiquke 14.— Residual stresses In turbine disk upon completion ol one operating cycle. 



APPENDIX 

STRESS CALCULATION FOR DISK WITH CENTRAL HOLE 


The calculations given in the section Numerical Examples 
deal with a disk that is solid at the center and has temperature 
gradients such that the plastic flow is confined to the region of 
the rim. Disks of other types spun under different conditions 
may be subject to plastic flow in other regions. 

One example of such a disk is a parallel-sided disk with a 



Figure 16.— Stresses In parallel-skied disk with central hole. 

central hole spun at a uniform temperature. In this disk 
plastic flow first occurs in the region of the central hole. Such 
a disk spun at a speed great enough to cause some flow near 
the hole is calculated here. 

The essential columns of the elastic calculation are given 
in table VIII (a) . Flow is indicated at stations a, 2,3, and 4. 
However, as flow occurs the stresses farther out in the disk 
may be increased. The quantities depending on disk dimen- 
sions, together with the first approximation, shown in table 


VTTT (b), are found in the manner given in the text. When 
the values of B r , n and B l-n (columns 57 and 58) are found for 
the stations at which flow occurs, however, new values of 
B T , n and j B 1>a must also be computed for all other point sta- 
tions before a new value of (column 59) can be found. 
These computations are also shown in table VILL (b) for the 
first approximation. Additional approximations must be 
made in the same manner until the correct flow increments 
are found. The stresses so calculated are plotted in figure 16. 

Where large numbers of computations involving plastic 
flow at the center of the disk are to be made, it may be desir- 
able to change the finite-difference approach to the problem 
in such a manner that the calculations are made from the 
outside of the disk toward the center instead of from the 
center toward the rim. This procedure has certain disad- 
vantages as a general approach to the problem of stresses in 
disks, particularly in that it requires a greater number of 
significant figures to obtain the same accuracy. For special 
applications it may, however, present- advantages that out- 
weigh the disadvantages in more general problems. 
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TABLE I.— CALCULATION OF STRESS DISTRIBUTION FOR EXAMPLE I 

(a) ElasttMtress calculation 
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TABLE I.— CALCULATION OF STRESS DISTRIBUTION FOR EXAMPLE I— Concluded 


(a) Elastic-stress calculation— Concluded 


SS 



0 

-15.292 
-38.834 
-95.481 
-172.39 
— 268.85 
— 509.83 
— 1, 287.3 
-3. 478. 6 
-4*913.8 
-8,997.6 
—11X469 
— 16,366 
— 26,915 
-46,028 
-79,644 
-106,070 
-139,980 
-181,250 
-213, 3S0 


37 


ri.„ 

[s-r.j-iseid 
+ ( 33 }i 


Constant at 
27,755 for 
ffr.SOf 4595 
Ib/sqln. 

(flg.fl 


(33)X(37)- 

( 36 ) 


27,766 
27,719 
27,578 
27,576 
27,447 
27,292 
28,893 
29,076 
31,608 
35,222 
37,779 
38,300 
36,689 
33,843 
23,598 
18,702 
16, 412 
13,374 
9,693 
4,595 


39 




27,756 
27,740 
27,718 
27,659 
27,682 
27,488 
27,246 
28,108 
29,669 
30, 932 
81,493 
29,952 
25,103 
15,131 
-6,604 
-41,628 
-67, 593 
-101,090 
-141,720 
-174, 760 


40 


73,500 

73,500 

73,500 

73,600 

73.600 
73,500 

78.600 

73.600 
73,600 
73,600 
73,500 
73,500 
73,600 
73,600 
73,500 
73,600 

72.000 
70,800 

67.000 

54.000 


41 


J CS 8 )H-( 39 )i 

y-mxm 


27,766 

27,730 

27,608 

27,618 

27,616 

27,391 

27,071 

28,604 

30,636 

38,286 

86,061 

84,883 

82,485 

29,364 

27,632 

63,490 

77,120 

108,400 

146,760 

177,090 


(b) Constants determined by disk geometry to be used 
In all plastic calculations 


n 

42 

43 

44 

46 

0»/(i ) 

CMOrt 

1-H42) 

i-c«) 

17. 

18. 
19- 
b„ 

a 016152 
.014706 
.014286 
.013889 

0.015626 

.016162 

.014706 

.014286 

L0152 

L0147 

L0143 

1.0139 

0. 9S43S 
. 96435 
.98529 
.98571 


(c) Plastic-stress calculation 


Approximation 

R 

48 

46a 

47 

48 

49 

60 

51 

52 

53 

*•>» 

(graph 
reading, 
fig. 7) 

(estimate) 

A*,*., 
(46a) 
(82)X2 A 
12(80) -(61)1 

Aj,a, 

(4«a) 

(82) X2 
[2(81) -(80)] 

(47) X(42) 

(47)»-iX(43) 

(48)X(44) 

(4S)--iX(45) 

Jpr 

(49)+(&)- 

(61>+(52> 

1 

17 

20X10-* 

20X10-* 

13.021X10-* 

-19.657X10-* 

Ol 19729X10-* 

0. 00000X10-* 

-19.956X10-* 


20. 153X10** 


18 

670 

670 

395.07 

-066.15 

5. 809C 

.19729 

-676. « 

-19.359 

662.69 


19 

I860 

1850 

1014.2 

-1817.1 

14 489 

6.8099 

-1873. 5 

— 658*35 

1237.4 


b 

3960 

3960 

2056.6 

-3969.1 

28.664 

1A 489 

-4014.1 


223 a 5 

2 

17 

20 

20 

12.582 

-19.769 

. 19034 

0 

-20.059 

o 

20.249 


18 

690 

6S5 

402. 14 

-881.31 

A 9139 

.19034 

-69^33 


677.97 


19 

1880 

1865 

1040.3 

-1SM.I 

14.862 

A 9139 

-188^3 

-071.29 

123A8 


b 

4650 

4300 

2312.1 

-4295.8 

32.113 

1A 862 

-4355.6 

-183A1 

2508.4 

3 

17 

20 

20 

12.544 

-19.782 

.19007 

0 


0 

20.252 


18 

670 

670 

393.05 

—686.43 

A 7802 

.19007 

-676.23 

-19.463 

662.74 


19 

1850 

I860 

1036.5 

-1855.8 

It 807 

5.7802 

-1882.3 

- 66 a 63 

1246.3 


b 

1900 

4100 

2218-3 

-4096.2 

30.810 

1A 807 

-4162.1 

-1829.3 

2368.4 

4 

17 

20 

20 

12.665 

-19.780 

.19023 

0 


0 

20.250 


18 

690 

680 

398.99 

-676. 38 

6.8675 

.19023 

-686.30 

-19. 461 

672.90 


19 

I960 

1880 

1037.2 

-1856.7 

14.817 

. 5. 8676 

-1882.2 

- 66 a 42 

1236.5 


b 

2900 

4000 

2158.0 

-3995.8 

29.946 

1A 817 


-1829.2 

226A9 

5 

17 

20 

20 

12.668 

-19.759 

.19028 

0 

—2a 059 

0 

20.249 


18 

680 

680 

399.27 

-678.33 

A 8717 

.19028 

-688.27 


672.87 


10 

I860 

1860 

1037.5 

-1865.7 

14.822 

A 8717 

-1882. 2 


1236.5 


6 

4200 

4015 

2160.3 

-4011. 0 

30.004 

14.822 

-4066.8 

— 1S29.2 

2282.4 

6 

17 

20 

20 

12.557 

-19.700 

.19026 

0 


0 

20.260 


18 

680 

6S0 

399.26 

-076.33 

5.8714 

.19026 

-686.27 

-19.461 

672.87 


19 

1860 

1860 

1037.6 

-1856.7 

14.822 

A 8714 

-1882.2 

-666.38 

ime 


5 

4015 

4015 

2160.8 

— 4011.0 

3 a on 

14.822 

— 106A8 

-1829.2 

2282.4 ... 


Approximatton 

ft 

64 

55 

66 

57 

58 

59 

60 

81 

62 

WSt 

(25b) -(53) 

AT.. 

(54) X (10) + 
(14)X(21) + 

AT., 

(54)X(8) + 
(14>X(20) + 
(29) 

Sixfsffclt 

(Si) 

n9)x©^+ 

(30)X«8)^i+ 

(56) 

<rt., 

6rr,t-(57)l] 
. +(33)» 

(WX(39>+ 

C34)X(69)+ 

m 

Vi.a, 



n 


-3174.2 

—37, 271 

-20,718 

-70,644 

27,755 

18,702 

-41,628 

63,490 

calculation, table I (a) . 

17 


-1538.9 

-28,938 

-25,257 

-10a070 



16, 412 

— fli, 593 

77, 120 


18 

1407.2 


-38,664 

-31, 137 

-139,880 


I A 374 


108,400 


19 



-49,881 

-39,667 

-181,250 


9,502.7 

-141,720 

14a760 


b 



-68,543 

-60,232 

-213,380 



4,595.2 

IrfMt i.U 

177,090 

1 

o 



-28,374 

-25,249 

-105,600 

28,023 

13,819 

-69,424 

77,280 



74A61 

-1444.6 


-30,848 

-121,480 

in 


-85, 019 

90,0® 


Era 


- 145 a 6 

-18,164 

-38,266 

-132,460 


7,825.9 

-9A423 

99, 507 


b 

427.90 

-137a 0 


-46,810 

-124,080 


4,695.8 

—87,881 

90,247 

2 

17 

999.45 


-28, 371 

-26,249 

-105,600 



— 69, 616 

77,295 



729.23 

—1438.5 


-30,842 

-12U070 


10,785 


90,577 


BrB 


-1456.8 

-18, 169 

-33)248 

-132,130 


7,727.0 

-95,147 

99,237 





—2 544.3 

—4a 680 

—116,510 


‘ 4,695.5 

-8a413 


3-_ 



-163a 3 

-28,371 

-25,219 

—105,600 

2A9S8 

13,782 

-69,459 

77,277 



744.46 

-1444.5 

-2a 698 

-80,348 

-m, 4 so 


10.846 


90,963 




-1453.3 

-17,925 

-38,263 

-132,260 


7,784.7 


99,349 




-1336.1 

-a 904. 9 

-46,761 

-121,010 



4,595.4 

—84,826 

87,216 

4 


999.45 


-28,371 

-25,249 

—105, 500 

26,008 

13,797 

—68,445 

77,273 






-30,844 

-121,200 


10,866 

-84,760 

90,682 



ML 90 

-145a 9 

-18, 177 

-38,254 

-182,260 


7,811.4 

-9A2Q4 

99,340 


HMH 

■MM 

-1366.8 

-0,117.9 

-4a 781 

-123,220 


4,695.2 


89,408 

6 


999.46 


-28,371 

—25,1149 

-10A500 

•28,006 

13,794 

-69,448 

77,274 



73A33 

■ICE™: 


—30,844 

-121,200 



-84,763 

SO, 683 



69L9G 

-145a9 

-18,177 

-38,254 

-132,250 



-95,207 

99,341 


b 


-1362.1 

-8)779.9 

-43,777 

-122,830 

— 

4,595.3 

-86,655 

89,072 


1 This value of <rt,« is also substituted for the original value of used In table 1(a) to compute plastic stress for stations a to 16. 
005385—50 20 
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TABLE II.— ABRIDGED VALUES FROM CALCULATION OF STRESS DISTRIBUTION FOR EXAMPLE II 



TABLE III.— ABRIDGED VALUES FROM CALCULATION OF STRESS DISTRIBUTION FOR EXAMPLE III 



TABLE IV.— CALCULATION OF STRESS DISTRIBUTION FOR EXAMPLE IV 

(a) Abridged values 
































































DETERMINATION OF STRESSES IN GAS-TURBINE DISKS SUBJECTED TO PLASTIC FLOW AND CREEP 


291 


TABLE IV.— CALCULATION OF STRESS DISTRIBUTION FOR EXAMPLE IV— Concluded 


(b) Calculation at effect of oreep on final stresses 


21 

63 

04 

66 

06 

67 

08 

68 

70 

71 

Cn- 

TV, 

(63) X5 

(64) v 
(40X2* 
(2X(3S) -036)1 

it.., 

(64) 

(41) X2* 
[2X (38)— (33)] 

(65)X(42) 

(63»+X(43) 

(66)X(44) 

(66)^.X(45) 

O'.. 

(67) +(68) — 
(68) +(70) 

17- 

18. 

19. 

6- 

0.2XKH 

.8 

.2 

0 

1.0X10"* 

1.5 

1.0 

0 

0.70226X1IH 

1.0848 

.80246 

0 

-0. 90760X10-* 
-1.4396 
-.81800 
0 

0. 010641X10-* 
.016953 
.011464 

0 

0.00000X10"* 

.010641 

.015653 

.011464 

— 0. 8S237X10 -1 
-1.4608 
-.83113 
0 

0.00000X10-* 
— .953C0 
-1.418-4 
-.60483 

0. 99301 X1(H 
.53439 
-.45935 
-.89342 

n 

n 

73 

74 

75 

70 

i ! 

“*■ 78 

78 

Il'.-IP'J-lQ'.l 
-Q .. 
(25b)— (71) 

if.. 

[(72mi«+at) x 

if'., 

I(7»X(8)+(14)X 
(20)j +(28) 

(28)X(7S).-i+ 

(73) 

A,. 

( 2 fl)xm).-i+ 

(30)X(7ffl.-i+ 

(74) 

[r r ,»— (75) &]+ 
(33) k 

C33)X(77)+(76) 

(34)x'(77)+(7n) 

17. 

18. 
10. 
&- 

638. 54X1CH 
-405. 17 
—647. 04 
-1400.5 

— 2460. 1 
-2247.2 
-2300.2 
-2230.7 

-12,873 

7,607.5 

11,751 

23,606 

1 1 1 1 
$$$& 
3SS3 

-106, 160 
—01,828 
—70,873 
-50,665 

46,005 

16,105 

13,199 

11.003 

9,334 

-45.443 
-32, 418 
-10,536 
10, 470 


TABLE V.— ABRIDGED VALUES FROM CALCULATION OF STRESS DISTRIBUTION FOR EXAMPLE V 


■ 

n 

2 

5 

6 


25 

25a 

25b 

38 

mm 

54 

60 

61 

| 

■ 

a 

D 


AT. 

ITn 

[P'J+EC'J 

sv-IP'.I-N'J 

D 

trt.n. 

HV-[PV)-[Q’.l 


ft.n 

PC 

0.5000 

4.3750 

2a 800X10* 

aS2SQX10-* 

331 




22,285 

21285 


21285 

22.286 

2„ 

.6260 

4.3750 

2a 800 

8.8300 

332 

9.49X10 -4 

0 X10-* 

9.49 X10-* 

22,253 

21023 

9. 49X10-* 



3.. 

.7800 

4.3780 

28.800 

a 8310 

333 

9.16 

0 

9.16 

22,181 

21,811 

9.10 

21181 

21,811 


■ It | <!(■ 

4.3750 

28.800 

a 8330 

334 

a 50 

0 

9.50 

21027 

21,645 

9.50 

21027 

21, 645 

*~ 

L2500 

4.3750 

28.700 

8.8390 

338 

37.36 

0 

37.36 

21,801 

20,887 

37.36 


20,687 

6.. 

LJ000 

4.3750 

28.700 

8.8440 

341 

2A22 

0 

2A22 

21,501 

20,085 

2a 22 

21,501 


■trl 


4 3750 

28.700 

8.8580 

350 

84.50 

0 

84.60 

2A 744 

1A 160 

84.50 

20,744 

gaBTCTgB 

B: 1 

Bin « B 

A8400 

28.500 

A 8970 

374 

226.70 

0 

22A70 

21,090 

14,149 

22a 70 


HUB 

B ■ 


3.2750 

28.300 

A 9530 

409 

334.78 

0 

334 78 

20,964 

A 300 

334 78 



B ji 

■TVmB 

2.6800 

28.100 

9.0240 


42A09 

0 

428.09 

20,680 

1,010 

42A09 

20,680' 

1, 010 

B[ j 

BjJ/rriB 

2.3720 

28.000 

B.06S3 

479 

254.27 

0 

254.27 

20,585 


254.27 

20,685 

-3,342 

Bn 

ByfflB 

2.2100 

27.800 

9.1120 


287.76 

0 

287.76 

19, 144 


287.76 

10,144 

-A 724 

13_ 

6.6000 

2.1600 

27.600 

0.1020 

539 

309.42 

0 

309.42 

16,529 

-14,887 

309.42 

16, 629 

— 14,887 

KS 


2.1550 

27.400 

0. 2140 

572 

332.09 

0 

'33109 

13,400 

-21,552 

332. 09 


-21,552 

Ki 

B*y»>B 

2.7000 

27.100 

9.2720 


366.97 

0 

30a 97 

7,748 

-29,867 

3C0.97 

7,748 

-29,867 

16. 

8.0000 

2.7000 

20.800 

9.3330 

646 

39L74 

0 

391.74 

4,365 

-37,781 

39L74 

4,365 

-37, 781 

17 

8.2500 

15500 

2a 700 

0. 3660 

066 

207.97 

107. 16 

40.840 

A 879 

-37,429 

40.840 

2,87# 

-37,429 

18. 

8.6000 

2.3800 

2a 500 

9.3900 

687 

220.02 

858.07 

-038.97 

1,076 

-19, 393 

-63a 07 

1,675 

-19, 393 

19. 

8.7500 

11450 

26.300 

9.4490 

70S 

232.78 

1338.5 

-1100. 5 

1, 171 


-1100.5 

1, 171 


6- 

9.0000 

L 9100 

20. 100 

9. 4990 

730 

244.38 

21015 

-1918.5 

1,704 

58,858 

-1918.6 

1,764 

51858 


TABLE VI— ABRIDGED VALUES FROM CALCULATION OF STRESS DISTRIBUTION FOR EXAMPLE VI 



i 

2 

5 

6 

1 

25 

25a 

25b 

38 

39 

54 

00 

61 

n 


B 

E. 

a*. 

AT, 

H’. 



«>.» 







4.3750 

29.300X10* 

a 6700X10-* 

230 




3616 

3,616 


3760 

a 750 

2.. 

.6250 

4 3760 

29.300 

A 6770 

237 

McrW'iUilii B 

0 X10-* 

A 91X10-* 

3592 

A 37* 

a 91X10-* 

3726 


3., 

.7500 

4.3750 

29.300 

a 6700 

238 

9.15 

0 

9.15 

3637 

A 160 

9.15 

3672 

3,294 

n 


4.3750 

29.300 

a 6840 

241 

27.20 

0 

27.20 

3364 

2,529 

27.20 


2,662 

6 .- 

L2500 

4.3750 

29.200 

a 6880 

244 

27.03 

0 

27.03 

3140 

1,956 


3274 


6.. 

LS000 

4.3750 

29.200 

a 6930 

247 

27.30 

0 

27.30 

2897 

1,402 

27.30 


1,536 

l 1 

mil 

4.3750 

29.200 

a 7010 

252 

4a 48 

0 

4a 48 

2417 

657 

4a 48 

2551 

GOO 

B; j I 

Bl nVjB 

3.8400 

29.100 

a 7190 

263 

100.44 

0 

100.44 

1725 

-1,301 

100.44 

1876 

-1,160 

B> i 1 

B2 

12750 

29.100 

A 7370 

274 

100.84 

0 

100.84 

933 

-A 144 

KHX84 

1100 

-2,990 

B • | 

Bl jwjB 

2.6800 

29.000 

a 7660 

2S6 

110.28 

0 

110.28 

2 

-5,260 

110.28 


-5,087 

B ' I 

B3 TjYjB 

13720 

29.000 

a 7640 

291 ' 

4a n 

0 

46.11 

-646 

-6,093 

46.11 

-317 

-5,907 

^p ' i 

Bi 53B 

12100 

28.900 

8.7730 

297 

5a 26 

0 

5a 26 

-1108 

-7,185 

6a 26 

-m 

-6,991 

B ■ 1 

m jjjB 

11600 

2a 900 

8.7820 

302 

4a 58 

0 

4a 68 

-1639 

-A 043 

4a 58 

-1393 

-7,843 

Brfl 


11550 

28.900 

a 7910 

308 

5a 46 

0 

6a 46 

-2140 

-9,150 

6a 40 


—8, 947 

15- 

7.6000 

17000 

28.900 

A 7990 

313 

40.40 

0 

4a 46 

-2165 

-A 819 

46.40 

-1974 


Kja 

b : 3v;vb 

17000 

2a 800 

a 8000 

319 

5a 98 

0 

55.98 

-2670 

—10,889 

sa 98 

-24S6 


Brfl 


15600 

2a 800 

a 8140 

322 

2a 04 

167.16 

-139.12 

-3025 

-0,756 

-139. 12 

-2824 

-a 604 

mm 

B:4'i"B 

13800 

28.800 

a 8170 

324 

18.60 

85a 97 

-810.37 

-3000 

17,132 

-84a 37 

-2785 

17,329 

mm 


11460 

28. S00 

8.8220 

327 

28.09 

1333.6 

-130a4 

-2197 

53,634 

-isoa4 

-1969 

53,740 

b.. 

9.0000 

L 9100 

2A800 

a 8260 

330 

27.79 

2162.5 

-2134.7 

0 

112,000 

—1461.0 

a 

9A673 
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TABLE VII— ABRIDGED VALUES FROM CALCULATION OF STRESS DISTRIBUTION FOR EXAMPLE VII 


It 

1 

2 

5 

6 

' 7 

25 

25a 

25b 

35 

39 

54 

60 

61 

r» 

h. 

E. 

a. 

A Tn 

H'„ 

[RJ+IQ'J 

H’.-tP'd-IC'J 


trt.n. 

ff.-tP'.j- 

fO'.l-PL 

<Tt,% 

<n.* 

a_ 

2.. 

3.. 
4__ 
6_. 

6.. 
7.. 
8_» 
9__ 
10. 
11- 
12. 

13. 

14. 

15. 

16. 

17. 

18. 
19. 
ft__ 

0. 5000 
.6250 
.7600 

1. JD000 
1. 2600 
1.6000 
2.0000 

3. 0000 

4.0000 

5.0000 
5.6000 

6.0000 

6.5000 

7.0000 

7.5000 

8.0000 
8.2500 
a 5000 
8. 7500 
9.0000 

4.3750 
4. 3750 
4.3750 
4.3750 
4.3750 
4.3760 
4.3780 
3.8400 
3.2750 
2.6800 
2.3720 
2.2100 
2. 1600 
2.1550 
2. 7000 
2.7000 
2.5500 
2.3800 
2.1460 
1.9100 

§ 

X 

s 

"S 

I 

Constant at 8.2970X10** 

© 

s 

1 

i 

D 

0 

■J 

1 

0 xio^ 
Q 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

167. 16 
858. 97 
1333.6 
1488.8 

0 X1Q-* 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

-167. 16 
-858.07 
-1333.5 
-1488. 8 

-2468 

-2468 

-2468 

-id® 
■-268 
-2J88 
—2468 
—2772 
-3183 
—3775 
-4202 
-4437 
— 4476 
-4432 
-3518 
-5615 
-$44 
-3B78 
-3282 
0 

-2,468 
-2,408 
-2,468 
-2,468 
-2, 468 
-2,468 
-2,468 
-2,614 
-2,849 
-8,193 
-8,428 
-3,004 
-3,705 
-3,762 
-3,476 
-3,479 
1,358 
26,512 
65,023 
107,180 

6 xich 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-167. 16 
-858.97 
-1333.5 
-947.81 

ii m i i i i i i i i i i m i i i 

IIIIIIHIIMIIM 


TABLE VIII.— CALCULATION OF STRESS DISTRIBUTION FOR PARALLEL-SIDED DISK 

(a) Abridged values 


n 

1 

2 

5 

6 

7 

25 

25a 

25b 

38 

39 

40 

41 

60 

01 

r m 

h. 

E, 

a» 

AT, 

H\ 

tP'J+iQ'J 

H’.-tP'.l-iQ'J 

*>,« 

*t.n 


<r«,» 

+r.» 

«*la» 


0.6000 








0 

176,040 


176,040 

0 

124,620 

2.. 

.0250 








31,650 

142,520 


129,620 

24,062 

127,160 

3.. 

.7500 








48; 400 

124,850 


109,030 

41,518 

124,510 

4- 

L0000 


%. 

t 





64,573 

100,350 


92, 802 

60,031 

111,630 

6 „ 

1.5000 

_ 


s 





74,394 

91,900 

X 

84,518 

72,549 

93,977 

6.. 

2.0000 

g 

X 

X 

o 

o 

O 

O 

75,607 

86,678 


81,704 

74,625 

87,813 

7.- 

2.5000 


§ 

p 


"S 

H 

05 

74,187 

83,268 


79,119 

73,589 

84,005 

8__ 

3.0000 



Q 

+J 

"3 

+» 

+3 

71,369 

80,278 

•5 

76,215 

70,976 

80,803 

9„ 

3.5000 

«s 


cc 

a 

S 

9 

9 

67,553 

77)268 


72,898 

67,282 

77,669 

10_ 

4.0000 

ri 


3 


a 

ta 

ta 

62,892 

74,000 

g 

69,156 

02,700 

74,380 

11. 

4.5000 

3 





O 

O 

57,471 

70,579 

S 

65,024 

67,332 

70,845 

12. 

5.0000 

a 


| 

o 

O 

O 

O 

51,322 

06,782 

3 

60,551 

51,223 

67,009 

IS. 

5.5000 

o 

g 






44,473 

62,650 


65,827 

44,401 

02,847 

14. 

6.0000 

O 

* g 

I 





36,934 

58,163 


50,979 

36,883 

8*330 

15. 

6.5000 


o 

D 





28,711 

53,317 


46,219 

28, G78 

53,476 

10. 

7.0000 








19,813 

48,103 


41,873 

19, 792 

48,249 

17. 

7.6000 








10,248 

42,518 


38,433 

Id 232 

42,654 

b.. 

8.0000 




. 




Q 

80,562 


36,592 

0 

36,91# 


(b) Calculation for first approximation of plastic-stress distribution 


n 

42 

43 

44 

45 

46 

4fla 

47 

48 

_ 40 

50 

61 

52 

53 

0»/<l) 


1+C42) 

1— (43) 

(grapji read- 
ing, flg. 7) 

S.« 

(estimate) 

Ar.», 

(46a) y 
1417X2* 
PX (38) -(39)1 

A*,*, 

(46a) " 
(SJ><2 X 
[2X(39)-(38)] 

■ (t7)X (42) 

(47) .-iX(43) 



(49)+JsO)- 

w+m 

2.. 

8.. 

4_. 

mgmm 


B jm 1 ! 

B4US1 

a 87500 
.90000 
.83333 

HFvjjH 


-900.00X10"* 

-168.00 

-38.59 

0 

1800.00X10^ 

587.56 

276.95 

0 



M ‘T r a ' rtaial 


■mppi 

HI 

-18. 806X10-« 
-A 2158 
, 0 










































